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A search for muon neutrinos from neutralino annihilations in the Sun has been performed with 
the IceCube 22-string neutrino detector using data collected in 104.3 days of live-time in 2007. 
No excess over the expected atmospheric background has been observed. Upper limits have been 
obtained on the annihilation rate of captured neutralinos in the Sun and converted to limits on the 
WIMP-proton cross-sections for WIMP masses in the range 250 - 5000 GeV. These results are the 
most stringent limits to date on neutralino annihilation in the Sun. 
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Non-baryonic cold dark matter in the form of weakly 
interacting massive particles (WIMPs) is one of the 
most promising solutions to the dark matter problem [l[ . 
The minimal supersymmetric extension of the Standard 
Model (MSSM) provides a natural WIMP candidate in 
the lightest neutralino Xi ■ This particle is weakly in- 
teracting only and, assuming R-parity conservation, is 
stable and can therefore survive today as a relic from the 
Big Bang. A wide range of neutralino masses, m^o, from 
46 GeV 0] to a few tens of TeV 0] is compatible with ob- 
servations and accelerator-based measurements. Within 
these bounds it is possible to construct models where the 
neutralino provides the needed relic dark matter density. 

Relic neutralinos in the galactic halo may become grav- 
itationally trapped in the Sun and accumulate in its 
center, where they can annihilate each other, produc- 
ing standard model particles. These may decay, cre- 
ating neutrinos which can escape and reach the Earth. 
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The search presented here aims at detecting neutralino 
annihilations indirectly by observing an excess of such 
high energy neutrinos from the Sun. Limits on the neu- 
trino flux from the Sun have previously been reported by 
BARS AN 0, MACRO 0, Super-Kamiokande 0], and 
AMANDA [|. 

The IceCube detector [§] records Cherenkov light in 
the ice from relativistic charged particles created in neu- 
trino interactions. In 2007 the detector consisted of an 
array of 22 vertical strings with 60 Digital Optical Mod- 
ules (DOMs) each, deployed in the clear Antarctic ice at 
the South Pole at depths between 1450 m and 2450 m be- 
low the ice surface. The vertical spacing between DOMs 
is 17 m and the horizontal distance between strings is 
125 m. Each DOM consists of a pressurized glass sphere 
containing a 25 cm photomultiplier tube (PMT) and a 
digitizer board. The PMT waveforms are stored when 
nearest or next-to-nearest DOMs fire within 1 /is. The 
trigger selects time windows when eight DOMs produce 
waveforms within 5 /is. The reconstructed first photon 
arrival times are used to determine the muon direction. 

The background in the search for neutrinos from the 
Sun comes from air showers created by cosmic ray in- 
teractions in the atmosphere. The showers cause down- 
wards going atmospheric muon events, triggering at sev- 
eral hundred Hz, and atmospheric muon neutrino events, 
triggering at a few mHz. When the Sun is below the 
horizon, the neutrino signal can be distinguished from 
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FIG. 1: The product Qi x Q2 of the output values of the two 
SVMs for the experimental data, a simulated signal (m^o = 
1000 GeV, hard spectrum) and the background. The back- 
ground has been scaled to match the data rate and it is shown 
divided into three components: atmospheric neutrinos and 
single and coincident atmospheric muons. 
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FIG. 2: Cosine of the angle between the reconstructed track 
and the direction of the Sun, for data (squares) with one 
standard deviation error bars, and the atmospheric back- 
ground expectation from atmospheric muons and neutrinos 



(dashed line). Also shown is a simulated signal (m^- 



GeV, hard spectrum) scaled to [i 3 



1000 



events (see Table I). 



the atmospheric muon background by selecting events 
with upward-going reconstructed muon tracks. 

The dataset used in this analysis consists of ~ 4.8 TO 9 
triggering events taken while the Sun was below the hori- 
zon, corresponding to 104.3 days of livetime between 
June 1st and September 23rd, 2007. The events were 
processed through several filters to reduce the content of 
atmospheric muon events and to enrich the dataset in 
muon- neutrino events. The analysis was performed in a 
blind manner such that the azimuth of the Sun was not 
looked at until the selection cuts were finalized. 

Events were first required to have at least ten hit 
DOMs, and the zenith angle of the line-fit [T3| first- 
guess reconstructed track was required to be larger than 
70°. Selected events were subjected to Log-Likelihood 
(LLH) fitting of muon tracks [lfj| . which uses the prob- 
ability distribution of the photon arrival times. Cuts 
were then placed on the zenith angle of this reconstruc- 
tion (90° < 6> LLH < 120°) and the width of the likeli- 
hood optimum (ctllh < 10°), to select upwards going 
events of good quality. Very loose cuts were placed on 
several kinematic quantities to remove a small number 
of outlying events. The final background reduction was 
then done using Support Vector Machines (SVMs) [ll[, 
multi-variate learning machines used to classify events as 
signal-like or background-like. Twelve event observables, 
that correlated modestly with one another (correlation 
coefficient |c| < 0.5), were used to train two SVMs with 
six input observables each. The use of two SVMs allowed 
minimal correlation (\c\ < 0.3) between the six observ- 
ables for each SVM. Training was done with simulated 
signal events, and a set of real data, not used in the anal- 
ysis, was taken as background. The observables describe 



the quality of the track reconstructions and the geometry 
and the time evolution of the hit pattern, most notably 
through the opening angle between the line-fit and the 
LLH tracks, ctllh, the mean minimal distance between 
the LLH track and the hit DOMs, and the number of hit 
strings. The SVM input distributions for data and sim- 
ulated backgrounds were generally in good agreement. 

Three types of background were simulated: atmo- 
spheric muon events from single and coincident air show- 
ers were simulated using C0RSIKA [T2], and atmospheric 
Ufj, events were simulated following the Bartol spec- 
trum [l3[ . Solar- WIMP signals were simulated with 
WimpSim (l4| . Two neutralino annihilation channels, 
W + W~ (hard channel) which produces a harder neu- 
trino energy spectrum, and bb (soft channel) which gives 
rise to a softer neutrino energy spectrum, were simulated 
for five masses m^o = 250, 500, 1000, 3000, and 5000 
GeV. The neutrinos were propagated through the Sun 
and to the Earth with full flavour oscillation. Absorp- 
tion in the Sun is important for neutrinos with energies 
above a few hundred GeV. A muon and a hadronic shower 
were generated in the ice near the detector. At the ver- 
tices the mean energy of simulated signal muons ranges 
from about 30 GeV to about 150 GeV depending on sig- 
nal model, see Table I For the hard channel < E t , > 
decreases for m^o > 3 TeV owing to neutrino absorp- 
tion in the Sun and secondary neutrino generation. The 
muon contribution from tau decay was evaluated to be 
insignificant and tau vertices were therefore neglected. 
Propagation of muons through the ice was simulated [l5| , 
and the Cherenkov light propa gati on from the muon to 
the DOMs was performed with fig , taking into account 
measured ice properties (l7| . 

Fig. [1] shows the distributions of the product of the 
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FIG. 3: Upper limits at the 90% confidence level on the 
muon flux from neutralino annihilations in the Sun for the 
soft (bb) and hard (W + W /_ ) annihilation channels, adjusted 
for systematic effects, as a function of neutralino mass. The 
shaded area represents MSSM models not disfavoured by di- 
rect searches [2ll . |22| . A muon energy threshold of 1 GeV 
was used when calculating the flux. Also shown are the limits 
from MACRO @, Super-K 0, and AMANDA [|. 



two SVM output values, Qi x Q 2 . As can be seen in 
the figure the distribution of simulated background is in 
good agreement with data. The final event sample was 
selected by requiring Qi x Q 2 > 0.1. This cut increased 
the signal : Vbackground ratio by a factor of 8. 

Simulations predict that the final data sample of 6946 
events has an atmospheric f M event content of 56%, and 
that the remainder consists of mis-reconstructed atmo- 
spheric muon events. The loose cuts maintain a large ef- 
fective volume, defined as the detector volume with 100% 
selection efficiency, since the final signal determination 
was done on the basis of direction. 

After calculating the Sun's position, the observed num- 
ber of events as a function of the angle to the Sun, ^, 
is compared to the atmospheric background expectation 
in Fig. [21 The angular distribution is consistent with the 
expected background and no excess of events from the 
Sun is observed. 

Using likelihood-ratio hypothesis tests the observed VP 
distribution is fitted with a sum of distributions of the 
simulated signal and the expected background. Here, 
the expected background is detemined by using real data 
with randomized azimuth direction of the Sun. We then 
follow the unified Feldman-Cousins approach [l8| to con- 
struct the confidence intervals on the number of signal 
events fi s . The upper 90% confidence limit ranges be- 
tween /z s = 6.4 and fi s — 8.5 events depending on signal 
case, see Tabic HI 

Simulation studies were used to estimate the system- 
atic uncertainty on the signal effective volume V c g . Un- 
certainties in the photon propagation in ice and absolute 
DOM efficiency dominate, contributing ±17% to ±24% 



FIG. 4: Upper limits at the 90% confidence level on the spin- 
dependent neutralino-proton cross-section a SD for the soft 
(bb) and hard annihilation channels, adjusted for 

systematic effects, as a function of neutralino mass. The 
shaded area represents MSSM models not disfavoured by di- 
rect searches l2ll . |22| based on a SI . Also shown are the limits 
from CDMS [p, COUPP [H, KIMS 0] and Super-K @]. 



depending on the signal model. The total systematic un- 
certainty on V e fi ranges from ±19% for the highest m^o 
to ±26% for the lowest m^o. Deviations in the event rate 
between data and background simulations are within the 
systematic uncertainty. These uncertainties are included 
in the results presented below. 

From the upper limits on /x s we calculate the limit on 
the neutrino to muon conversion rate r = tt^t-, for 

" I* | — .Vrff-t ' 

the livetime t. Using the signal simulation [14|, we can 
convert this rate to a limit on the neutralino annihilation 
rate in the Sun, Ta, see Table Q] Results from different 
experiments are commonly compared by calculating the 
limit on the muon flux above 1 GeV, < i ) AI , which is also 
shown in Table U together with the sensitivity, the 
median limit obtained from simulations with no signal. 
A downward fluctuation in the data close to the position 
of the Sun results in limits lower than the sensitivity. 
Within "J < 3°, corresponding to the rightmost bin in 
Fig. [21 the fluctuation has a probability of 8.8%. In this 
bin we expect less than 0.4 background events from solar 
atmospheric neutrinos (l9| . 

The 90% confidence upper limit on $ M as a function 
of m^o is shown in Fig. [31 compared to other limits 
[10, Q, and MSSM model predictions 0. In the plot, 
the shaded area represents neutralino models not dis- 
favoured by the direct detection experiments CDMS (2lT | 
and XENON-10 [22[, based on their limit on the spin- 
independent neutralino-proton cross-section. 

The limits on the annihilation rate can be con- 
verted into limits on the spin-dependent, a SD , and spin- 
independent, a SI , neutralino-proton cross-sections, al- 
lowing a more direct comparison with the results of di- 
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TABLE I: Upper limits on the number of signal events fi s , the conversion rate F^^^, the neutralino annihilation rate in the Sun 
Fa, the muon flux $ M , and the neutralino-proton scattering cross-sections (spin-independent, a SI , and spin-dependent, a SD ), 
at the 90% confidence level including systematic errors. The sensitivity 3> M (see text) is shown for comparison. Also shown is 
the median angular error O, the mean muon energy <_B M >, the effective volume V e fj, and the effective area ^eff- 
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rcct search experiments. Since capture in the Sun is 
dominated by a SD , indirect searches are expected to be 
competitive in setting limits on this quantity. Assuming 
equilibrium between the capture and annihilation rates 
in the Sun, the annihilation rate is directly proportional 
to the cross-section. A limit on a SD is found by set- 
ting a SI to zero, and vice versa. We have used Dark- 
SUSY [2(| and the method described in [23| to perform 
the conversion. The results are shown in Table HI We 
assumed a local WIMP density of 0.3 GeV/cm 3 , and a 
Maxwellian WIMP velocity distribution with a dispersion 
of 270 km/s. Planetary effects on the capture were ne- 
glected. Fig. [4] shows the IceCube-22 limits on a SD com- 
pared with other bounds 0, HH, [H, [1| , and the MSSM 
model space defined as for Fig. [3] Indirect searches for 
dark matter in the Sun complement direct searches on 
Earth in several respects. WIMPs in the Sun would ac- 
cumulate over a long period and therefore sample over 
different dark matter densities in the galactic halo. This 
gravitational accumulation is sensitive to low WIMP ve- 
locities while direct detection recoil experiments are more 
sensitive at higher velocities. 



In conclusion, we have presented the most stringent 
limits to date on neutralino annihilations in the Sun, im- 
proving on the 2001 AMANDA [H limits by at least a 
factor of six for hard channels. We also present the most 
stringent limits on the spin-dependent WIMP-proton 
cross-section for neutralino masses above 250 GeV. The 
full IceCube detector with the DeepCore extension [26| is 
expected to test viable MSSM models down to 50 GeV. 
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